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Abstract 

One monolayer of Cgo on one monolayer of hexagonal boron nitride on nickel is investigated by 
photoemission. Between 150 and 250 K the work function decreases and the binding energy of the 
highest occupied molecular orbital (HOMO) increases by ~100 meV. In parallel, the occupancy 
of the, in the cold state almost empty, lowest unoccupied molecular orbital (LUMO) changes by 
0.4 lb 0.1 electrons. This charge redistribution is triggered by onset of molecular rocking motion, i.e. 
by orientation dependent tunneling between the LUMO of Cgo and the substrate. The magnitude of 
the charge transfer is large and cannot be explained within a single particle picture. It is proposed 
to involve electron-phonon coupling where C^q polaron formation leads to electron self-trapping. 
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In the emerging field of molecular electronics Ceo plays a pivotal role for the exploration 
and the test of new concepts. The separation of energy scales in inter- and intramolecular 
bonding gives rise to structural phase transitions that leave Cgo intact. This is e.g. seen 
with diffraction experiments where orientational ordering of Cqq was inferred from changes 
of the unit cell size Furthermore, Cgo compounds have narrow electron bands 

^, large on-site Coulomb repulsion and intramolecular phonons whose frequencies compare 
with the intermolecular electron hopping rates. The observation of high Tc superconductivity 
in alkali metal doped fuUerenes jsl or antiferromagnetism jgl emphasizes strong correlation 
effects. On this background a rich variety of thermodynamic phases emerges where in view 
of applications metal insulator transitions are particularly interesting. 

Here we report on the change of charge state of Cqq that occurs in parallel to an orien- 
tational order-disorder transition. The underlying substrate of the Ceo monolayer plays a 
crucial role: It provides the charge by tunneling, and it does not bind the Cqo molecules 
strongly. /i-BN/Ni(lll) is an ultimately thin, atomically flat, metal-insulator interface. Its 
geometry and electronic structure is well known from experiments [3, 1^ and density func- 
tional theory j^. Compared to Ceo/nietal interfaces 13, [l^, the substrate-induced doping, 
i.e. charge transfer on Ceo, is much lower which allows to investigate the insulator limit. For 
the description of the C6o//i-BN/Ni(lll) junction we borrow concepts from non-adiabatic 



gas surface reactions 



12| . Charge transfer onto a molecule occurs when the vertical affinity 



level is degenerate with the Fermi level of the substrate. Provided that the tunneling elec- 
tron resides long enough on the molecule, the molecule evolves towards its new equilibrium 
coordinates, i.e. approaches the adiabatic affinity level. The electronic system has to pay the 
energy for this Jahn- Teller distortion and can therefore not be described in a one-electron 
picture. 

The ex per iments were performed in a modified VG ESCALAB 220 photoelectron spec- 
trometer with an attached scanning tunneling microscopy (STM) chamber 
/i-BN/Ni(lll) interfaces were produced using the recipe by Nagashima et al. 
monolayers were prepared by annealing multilayers. On a linear heating ramp (0.15 K/s) 
Cqo desorbs from the multilayer at 475K and from the monolayer at 520K. From this it is 
inferred that Cqo is weakly bound on /i-BN/Ni(lll) 1.5 eV) [15]. At room temperature 
Cqo wets the substrate and forms a commensurate hexagonally close packed (4 x 4) structure 
with a Ceo nearest neighbor distance of 10.0 A (see Figure [T]). The STM picture in Figure 
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2Jd) shows how perfect the Cqq //i-BN/Ni(lll) junction grows. 

Figures ^ c) and d) are low energy electron diffraction (LEED) data [E = 23.5eV) at 
250 and 160 K, respectively. From comparison of the two it is seen that the unit cell size 
increases from a (4x4) with one Cgo molecule to a quasi (4V3x4V3R30°) unit cell with 3 Cqq 
molecules when the sample is cooled to 160K. This is in line with the findings of Benning 
et al. P and Goldoni et al. P], who found for Cgo multilayers in the same temperature 
range a (1x1) to (2x2) surface lattice transformation. This phase transition corresponds to 
an order disorder transition which occurs in bulk Cqq at :^250K ^|. 

In Figure El He la normal emission photoemission spectra for 300 and 160 K are shown. 
Due to the small mean free path of the photoelectrons the molecular orbitals of Cqq dominate 
the spectra and the underlying substrate contributes only weakly. The bump at 0.6 eV 
binding energy reflects the Ni d-band. The Cqq derived features are strongly temperature 
dependent. The highest occupied molecular orbital (HOMO) shifts from 2.75 up to 2.64 eV 
in going from 250 K to 160 K. Furthermore, the orbital peaks sharpen and the HOMO 
increases in intensity. From the large HOMO binding energy and the fact that there is 
no strong LUMO contribution which would produce a peak at the Fermi level we see that 
this system is an example for poor screening, i.e. high on-site Coulomb repulsion or large 
HOMO-LUMO separation in photoemission The work function increases parallel to 

the HOMO-shift from 4.25 to 4.37 eV. This indicates a charge redistribution in the interface 
and that the molecular orbitals are mainly bound to the vacuum level. 

The inset of Figure El shows the spectra in the vicinity of the Fermi level. Intriguingly, a 
80 meV down shift of the leading edge of the photoemission is observed in going from room 
temperature to 160 K. Such a shift is beyond instrumental uncertainties since the energy 
resolution was better than 60 meV. A careful comparison with the Fermi edge as measured 
on a Ag polycrystal indicates that the room temperature leading edge hes above Ep. As we 
show below it reflects a shift of the high binding energy wing of the LUMO parallel to the 
HOMO. The difference between the room temperature and the 160 K spectrum in the inset 
of Figure El corresponds to 4.5 % of the HOMO intensity at 300 K. This intensity can be 
assigned to the charge difference on the LUMO at high and low temperature, respectively. 

To get a more quantitative measure of the LUMO occupancy a one parameter fit was 
applied to the data. The fit procedure asserts the same Gaussian spectral shape as that of the 
HOMO and a degeneracy-derived weight of 10 for the hu HOMO and 6 for a fully occupied tiu 
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LUMO . The Gaussian for the LUMO is fitted to the spectrum that was divided by the 
corresponding experimental Fermi function (see inset of Figure Ej). The only fit parameter is 
the HOMO-LUMO energy difference , and the fit is performed in a region from the Fermi 
energy up to 200 meV above Ep. In order to obtain the LUMO occupancy the Gaussian for 
the LUMO is multiplied with the Fermi function and integrated. Figure El shows this LUMO 
occupancy versus the reciprocal temperature. It increases by a factor of 7 ± 1 in going from 
150 to 300 K Though these data are recorded during cool down of the sample, the 

effect is reversible with a hysteresis < 50 K. From the fit the LUMO occupancy changes 
between room temperature and 160 K by ~ half an electron. This is in good agreement 
with 0.45 electrons [2M indepenently derived from the difference of the spectra near the 
Fermi energy (inset of Figure |2I) and the comparison with the HOMO intensity. In Figure 
El we show as well the intensity of the HOMO. It decreases in a non Debye- Waller fashion 
with increasing temperature and indicates changes in molecular orientation as reflected in 
ultraviolet photoelectron diffraction (UPD) [2^. The transition temperature coincides with 
that of the change in charge state. This is a clear indication that the charge state is closely 
related to the molecular orientation. Since the spectra were recorded in normal emission it 
can be inferred that the change in orientation must contain rocking motion, i.e. rotational 
axes that do not coincide with the surface normal. 

How can we explain the change in charge transfer with temperature? In Bardeen's tun- 
neling picture the charge transfer between the substrate and the adsorbate is related to the 
overlap of the wave functions at the Fermi level [^|. The charge on the Ceo molecules thus 
reflects a finite tunneling probability of conduction electrons across the /i-BN layer onto the 
LUMO of Cqq. The residual occupancy of 0.08(2) electrons at low temperatures indicates 
that electron tunneling from the substrate to the molecule is still enabled, though at a low 
rate. At temperatures above 250 K the charge transfer begins to saturate. From this non Ar- 
rhenius behavior a model where the charge transfer is induced by thermally excited electrons 
alone, has to be refuted. However, the increase of the LUMO occupancy with temperature 
is related to the onset of molecular rocking motion in the Ceo monolayer (see Figure Ej) 
and therefore a molecular orientation dependent wave function overlap with the substrate 
gives a natural explanation for the effect. The onset of molecular rocking motion is in line 
with the known disorder transition in Ceo multilayers E], the LEED patterns in Figure 
X-ray photoelectron diffraction (XPD) and Ultraviolet Photoelectron Diffraction (UPD) 
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(also see Figure El) patterns. Furthermore, the XPD patterns (not shown) indicate that the 
Ceo molecules do not bind via the pentagons but rather via hexagons to the substrate, as it 
is e.g. the case for Cgo on Cu(lll) j22|. 

In order to justify the molecular rocking driven model, the electronic structure of Cgo 
has to deviate from spherical symmetry. For this purpose density functional calculations 
were performed within the Kohn-Sham scheme. As the exchange correlation functional we 
emp loyed the Generalized Gradient Approximation (GGA) of Perdew, Burke and Ernzerhof 
231 1. The Kohn-Sham orbitals were expanded either in a Gaussian basis set TZVP in the 

n 

TURBOMOLE code or in plane waves with a cutoff energy of 40 Ry. In the latter case 
the action of the core electrons on the valence electrons was replaced with norm- conserving 
pseudo potentials, and the length of the supercell was 21.2 A. 

Figure m shows results of density functional calculations of Ceo- In Figure E^) the degen- 
erate LUMO is shown as a constant electron density isosurface. Clearly, the LUMO electron 
density is largest on the pentagons, while the hexagons have a low density. Together with the 
fact that at low temperatures Ceo does not expose pentagons towards the substrate this gives 
a consistent picture: The overlap between the LUMO and the substrate, i.e. the tunneling 
rate increases if the molecules rock away from the equilibrium position and correspondingly 
more electrons tunnel into the LUMO. This overlap is shown in Figure |3Jd) as a function 
of the molecular orientation relative to the surface normal. The grey scale plot gives the 
integral of the LUMO in the half space z > Zo with Zq = 5A away from the center of Cgo 
as a function of the molecular orientation relative to this plane. The average probability to 
find an electron outside a plane z > Zo i.e. the mean of the integrals over all orientations 
/ \^LUMo{.z > -2o)p is 0.007 / I'^LUMol'^ cLud the anisotropy that we call tunneling anisotropy 
may be seen in the bottom panel of Figure |3)d). The tunneling anisotropy slightly depends 
on Zo, though for reasonable z^, > 5 A it never reaches values that could explain quantita- 
tively the strong temperature dependence of the LUMO occupancy. From the experiment 
the LUMO occupancy ratio between room temperature and low temperature is ~ 7, while 
a tunneling anisotropy of 0.75 could explain a ratio of 1.3 only. In other terms this theory 
explains an occupancy change of 0.01 electrons only, which is inconsistent with observation. 

Thus we have to refine the orientation dependent tunneling picture. It is unlikely that the 
substrate wave function and/or distance induced changes in the tunneling matrix element 
account for the strong amplification of the effect. Also, electron hopping between the Ceo 
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molecules is not expected to change the net charge balance in such a dramatic way. It is 
rather an indication that the one-electron picture as it is drawn from the tunneling anisotropy 
in Figure |3Jd) breaks down. The low tunneling rate provides substantial time for an electron 
to stay and to act on the Cgo layer before it returns to the substrate. For large residence 
times, comparable with the molecular vibration periods, the Cqq species start to deform 
and the electron couples strongly to phonons. In other terms, C^q distorts from the vertical 
affinity towards the adiabatic affinity E^. Our calculations on gas-phase Cq(^ result in a 
Jahn- Teller distortion energy, i.e. energy difference between unrelaxed and relaxed E^ — 
E^ of 51 meV. In the condensed phase we have to consider intramolecular ami extramolecular 
vibrations. Though, their interplay and relative importance is not known at the momen;. 
Intramolecular electron phonon coupling is e.g. observed in photoemission from {25 1, 
or recently in inelastic tunneling experiments on Cfin/AgfllO) , w here an electron energy 
loss of 54 meV was assigned to the excitation of a Hg phonon j26|. Near metallic surfaces 
the coupling to extramolecular vibrations occurs via the image force on the Cqq ion, which 
excites frustrated translational vibrations. In this image field the occupied LUMO dives 
below the Fermi level and gets trapped since back-tunneling to the substrate is prevented 
due to the lack of empty states. Together with the onset of Coulomb repulsion this may 
establish a new equilibrium LUMO occupancy. If indeed such a kind of 'polaron' driven self 
trapping mechanism is responsible for the strong amplification of the tunneling anisotropy 
this is a particular example for strong electron phonon coupling in such an interface system. 

In summary we report the observation of a temperature driven change in Ceo molecular 
orbital binding energies parallel to the work function and a strong change in LUMO oc- 
cupancy of Cqq on /i-BN/Ni(lll) that coincide with the onset of molecular rocking. The 
charge transfer is triggered by an anisotropic tunneling matrix element of the Ceo LUMO 
with the underlying substrate. The magnitude of the charge transfer is an indication for 
strong electron-phonon interactions. This effect might open the door for the application of 
the orientation dependence of the tunneling across molecules for electronic and spintronic 
switches. 
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FIG. 1: The investigated system C6o//i-BN/Ni(lll). a) Sketch of the interface, b) STM picture 
of the well ordered monolayer of Ceo- The three grey levels indicate three different terraces on the 
substrate. c,d) LEED patterns {E = 23.5eV): (4x4) at 250K and quasi (4^3x4^/3) R30° at 160K. 




FIG. 2: He la photoemission spectra of one monolayer of Cgo on /i-BN/Ni(lll) for room temper- 
ature and 160 K. The inset shows the strong shift of the photoemission leading edge that is due to 
the population of the LUMO at high temperatures. 
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FIG. 3: Occupancy of the LUMO (solid circles) and normal emission HOMO intensity (open 
triangles) as a function of reciprocal temperature. The solid lines are guides to the eye. The inset 
shows the extrapolated LUMO from the 300 K data in Figure [21 
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FIG. 4: Deviation of the LUMO of Ceo from spherical symmetry, a) Constant electron den- 
sity isosurface. b) Stereographic projection of the integral electron density outside planes placed 
5 A away from the center of Cgo- The low density parallel to hexagons (center of the plot) and 
the high density parallel to pentagons is apparent. The bottom panel shows the anisotropy of 
this electron density that is related to the tunneling rate between the substrate and the molecule 
along the dashed 9 line. The tunneling anisotropy of 1 corresponds to the average of all molecular 
orientations. 
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